Reentrant Phenomenon in Quantum Phase Diagram of Optical Boson Lattice 
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We calculate the location of the quantum phase transitions of a bose gas trapped in an optical 
lattice as a function of effective scattering length a e e and temperature T. Knowledge of recent 
high-loop results on the shift of the critical temperature at weak couplings is used to locate a nose 
in the phase diagram above the free Bose- Einstein critical temperature Tc°\ thus predicting the 
existence of a reentrant transition above T c , where a condensate should form when increasing a e g. 
At zero temperature, the transition to the normal phase produces the experimentally observed Mott 
insulator. 
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Optical lattices offer the intriguing possibility of inves- 
tigating properties of Bose-Einstein condensates (BECs) 
at varying effective interaction strengths Q, Q • Due to 
Bloch's theorem, particles loaded into such lattices be- 
have in some ways like particles in an ordinary continu- 
ous space. Thus they may form a BEC condensate in the 
zero-momentum state in the same way as in the contin- 
uum. 

Because of the lattice structure, the bosons in an op- 
tical lattice exhibit another phase transition which does 
not exist in the continuum [3|. It produces a state in 
which the atoms are frozen to their individual optical 
potential wells on the lattice and require an activation 
energy to be moved. This state is referred to as a Mott 
insulator due to its analogy to a state described long ago 
by Mott for fermionic systems For bosons on a lat- 
tice, a Mott insulator is generally expected to exist for 
integer filling factors at all temperatures if the repulsion 
between the atoms is larger than fegT so that the atoms 
have no place to escape from their optical potential wells. 

At very low temperature, the existence of this state has 
recently been demonstrated by determining the critical 
effective scattering length It is generally expected 
that at zero temperature the Mott transition should co- 
incide with the superfluid-normal transition. At the crit- 
ical point, the excitation energies of the bosons acquire 
a gap which pins the atoms to their potential wells. Ex- 
pressed differently, the Goldstone modes of translations 
have become massive and the associated phase fluctua- 
tions decoherent, in accordance with the criterion found 
in Ref. 0. 

Standard treatments of this situation proceed by 
analytic or numeric studies of the Bose-Hubbard model 
[B S The purpose of this paper is to simplify 
the model by focusing the attention upon systems with 
low-density where physics takes place in the lowest 
energy-band. There simple field-theoretic techniques 
can be applied to calculate the whole quantum phase 
diagram for the superfluid-normal transition. Our main 
point will be to include recent high-loop results near T c 
of the weakly coupled system to predict an interesting 
characteristic reentrant transition. In addition we shall 
observe that our calculated superfluid-normal transition 



point at zero temperature agrees with the experimentally 
observed Mott transition point. 

1) If bosons of mass M are trapped in a three- 
dimensional cubic periodic potential V(x) of lattice vec- 
tors S, i.e., F(x) = Vq X)j=i sin 2 ((7iXi) with qi = n/5, the 
wave vector q defines an energy scale E r = h 2 q 2 /2M. 
If the individual potential wells are deep, i.e., Vq ^> 
E r , the single particle Wannier functions u>(x) in the 
nearly harmonic wells are given by oscillator ground- 
state wave functions at the lattice sites 6 with size 
Aq = yJh/MujQ and energy hu « 2E r (V /E r ) 1/2 . The 
lowest energy band arising due to Bloch's theorem reads 
es(k) = 2J^ i=1 [l ~ cos(fci<5)] up to a trivial additive 
constant Here J follows from the tight-binding approxi- 
mation as J = / d 3 xw(x.)[-Ti 2 V 2 /2M + V(x)]w(x + 6) 

and is equal to J = iE r {V / E r f/ A e-^ v °/ E ^ n / ^ 
[9J- Due to the low-density of the system, a repulsive 
potential V(x) = gS^(x) with the coupling constant 
g = 4:TTh 2 a/M approximates well all relevant spherically 
symmetric short-range two-particle interactions, where a 
is the s-wave scattering length. In an optical lattice, this 
gives rise to an effective repulsive (5-function interaction 



with strength U = g J d 3 xw 4 (x) = {a / A$)2Tiw / \f2it — 

(27ra/A)VW^r(Vo/£r) 3/4 0, % The importance of 
the interactions between the particles in the periodic 
traps is measured by the ratio 7 = U/J between inter- 
action energy U = g e sn with g e ff = 47r7i 2 a e g / M e ff and 
kinetic energy J = h 2 n 2 / 3 /2M c g, where n is the parti- 
cle density (= f/S 3 with filling factor /). This leads to 
7 = 87ra G ff n 1 / 3 . 

The experimental optical lattice of Ref. is made 
of laser beams with wavelength A = 2<5 — 852 nm and 
contains about 2 x 10 5 atoms 87 Rb with a 1=3 4.76 nm 
[Tfl |. Its energy scale is E r w Ti x 20kHz «fcs x 150 
nK and Vo/E r is raised from 12 to 22. In this range, 
J/E r drops from 0.014 to 0.002, U/E r increases from 
0.36 to 0.57, and HoJo/E r increases from 0.36 to 0.57. Ex- 
panding the small-k behavior of the band energy £s(k) 
as ?i 2 k 2 /2M ff + . . . , the band width 4J defines an ef- 
fective mass M e ff of the particles M e g — h 2 /2JS 2 . In 
a typical BEC with a e e of the order of A and particle 
distances of a few thousands A, this ratio is extremely 
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small. For the particles tightly bound in an optical lat- 
tice, however, a e g n 1 / 3 can be made quite large. In the 
experiment of Ref.0, for temperatures near zero we have 
7 w 0.0248 exp(2 v / Vo/i?r), so that the increase of the po- 
tential depth Vo/E r from 12 to 22 raises a^n 1 / 3 from 1 
to 11.7. 

For increasing temperatures, we expect the critical 
a e ffn 1//3 to decrease until it hits zero as T reaches 
roughly the free BEC critical temperature Tj -* = 
27rft 2 [n/C(3/2)] 2 / 3 /M eff fc B with ((3/2) » 2.6124. In the 
above experiment where Vo/E r is raised from 12 to 22, 
the temperature drops from 14.2 nK to 1.93 nK, 
implying that T c (0) / ' E r drops from 0.094 to 0.013. Hence 
J and ksT are much smaller than Tiujq, so that we can 
ignore all higher bands. 

The purpose of this note is to derive the full tempera- 
ture dependence of this transition, thereby predicting a 
surprising reentrant phenomenon. 

2) We begin by considering a D-dimensional Bose gas 
in the dilute limit where the two-particle (^-function in- 
teraction is dominant. In the grand-canonical ensemble 
it is described by the Euclidean action 

AW,il>\ = dr d D xU*{x,T)[hd T + e(-ihV)-n] 

x^(x,t) + ^V(x,t)V(x,t) 2 }, (1) 

where \x is the chemical potential, and (3 = 1 /ksT . To de- 
scribe the phase transitions in this gas we calculate its ef- 
fective energy. We expand the Bose field T ) around a 
constant background i.e., ip(x.,T) = W + 5tp(x, r), and 
perform the functional integral for the grand-canonical 
partition function including only the harmonic fluctua- 
tions in Sip (x, t). This yields the one- loop approximation 
to the effective potential 
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with E(k) = V{e(k) -» + 2g cS |*| 2 } 2 - <? e 2 ff |*| 4 denot- 
ing the quasiparticle energies. An expansion parameter 
T] = 1 has been introduced whose power serves to count 
the loop order. The effective potential is extremized 
with respect to the background field 'J, yielding the con- 
densate density 

/i T) 2e(k) + [i 
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and the grand-canonical potential 
fi(/x,T) _ /i 2 
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The chemical potential fi is fixed by the total particle 
density n(/z, T) = — V~ 1 dfl(fi,T)/du. Eliminating n in 
favor of the condensate density no via 10, we find for the 
particle density 

e(k) + g eS n 



n-n = 
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= }+0(r, 2 ). (5) 

e /3Ve(k) 2 +2s„ ff n «:(k) _ 1 J 

This result is the so-called Popov approximation [Til IT^ . 
It is derived only for a small right-hand side where n « 
no. A standard way to extend such a relation to n > no 
is by making the equation self-consistent, replacing no 
by n (and rj by 1) on the right-hand side. The location 
of the quantum phase transition for all T is obtained by 
solving this equation for no = 0, 0] . The evaluation 
will be discussed in the next two paragraphs. 

Note that a more systematic approach to derive the 
self-consistent Popov approximation Q proceeds by 
applying variational perturbation theory to and 
(QJ according to the rules developed in I13 L | l 4l. and 
applied successfully to critical phenomena in 1 1 as well 
as many other strong-coupling problems |l6j |. There 
one introduces a dummy variational parameter jl by 
replacing ij, — > fi + rjr with r = (fi — p)/r], and re-expands 
consistently at fixed r up to the first power in rj. After 
this one re-inserts r = (n — jl)/rj and extremizes the 
resulting expression with respect to ji. The result turns 
out to be precisely the self-consistent version of the 
Popov approximation © (see also the related discussion 
of the self-consistent Hartree-Fock-Bogoliubov-Popov 
approximation in one-dimensional optical lattices 17J,|8|]). 

3) We first discuss the formation of a condensate for 
the free-particle spectrum e(k) = ?i 2 k 2 /2M G ff, where mo- 
mentum sums reduce to — > V J d D k / (2n) D . The in- 
tegral of the zero-temperature contribution can now be 
evaluated analytically, and we obtain in D = 3 dimen- 
sions for the transition curve in the T — a e g plane: 
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Here 1(a) abbreviates the integral 

f-OC> 

1(a) = I dx ■ 
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2^/xa + 16 (e v^W 1 ^ „ X ) 

and a = t 2 /a 2 s n 2 / 3 £(3/2) 4 / 3 is a dimensionless parame- 
ter, with t = T/Tc°^ being the reduced temperature. The 
result is shown in Fig. ^a). For small temperatures, the 
transition curve behaves like 

a cS n 1/3 = a + a x a + a 2 a 2 + O (a 3 ) , (8) 

with the dimensionless expansion coefficients a = 
(97T/64) 1 / 3 w 0.762, a x = -7r 2 a /24 w -0.3132, a 2 « 
0.1996, and w —0.207. The interaction causes an 
upward shift of the critical temperature from tc = 1 
to t c = 1 + 4 v / 2^v / aefin 1 / 3 /3C(3/2) 2 / 3 + O (a cS n^ 3 ). 
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FIG. 1: Quantum phase diagram of a homogeneous dilute 
Bose gas in variationally improved one-loop approximation 
without (dotted) and with properly imposed higher-loop slope 
properties at T c (dashed length increasing with VPT order 
from 1 to 3). Short solid curve starting at T = T c is due 
to a second-order finite-temperature calculation of Arnold et 
al. [Til . Dashed straight line indicates the slope of our curve 
extracted either from Monte-Carlo data |2(1 l2l| or precise 
high-temperature calculations 122.12311 . Diamonds correspond 
to the Monte-Carlo data of Ref. [24| and dots stem from Ref. 
|2*fij . both scaled to their respective values a e fs(T = 0). 



This has the square-root behavior found before in Refs. 

nana 

As announced in the abstract, the phase diagram has 
the interesting property that there exists a reentrant tran- 
sition above the critical temperature of the free sys- 
tem, which shows up as a nose in the transition curve, 
where a condensate can be produced by increasing a e s, 
which disappears upon a further increase of a e g. Our 
curves agree qualitatively with early Monte-Carlo simu- 
lations [24| as shown in Fig. ^ 

Recent Monte-Carlo simulation s Efl l2ll | and precise 
high-temperature calculations [22ll23|] indicate, however, 
that the square-root approximation is unreliable near 
Tc°\ the leading critical temperature shift being linear 
in the scattering length a c ff with a coefficient ci ss 1.3: 

t c = 1 + crflesn 1 / 3 + O (fleff" 2/3 ) • (9) 

It is possible to improve our self-consistent approxima- 
tion JHJ to accommodate the high- loop result 10. This 
can be done with the help of variational perturbation 
theory [H 03 El E| 

For this we use the expansion 
(JHJ) with a few exact coefficients and add two more trial 
coefficients to enforce the behavior 0. This produces a 
sequence of improved transition curves shown in Fig. ^ 
a) as dashed curves. 

4) We now consider an optical lattice where the 
wave vectors k of the band spectrum es(k) are re- 
stricted to the Brillouin zone ki S (— 7r/<5, 7r/5), and 
momentum sums become V YliLi J-i^/s dki/2n. The in- 
tegral is evaluated using the hopping expansion |29j . 
in which one expands the integrand in powers of the 
cosines in es(k). By doing so, we express the result in 



FIG. 2: Quantum phase diagram of optical boson lattice 
calculated for increasing hopping order (right to left). The 
quantity 7 = U / J is proportional to the effective scattering 
length a e ff. At zero temperature, the transition takes the 
superfluid into a Mott insulator. Near T c , higher corrections 
are needed due to infrared singularities. Since these are the 
same as in the continuum model, the correct curve should 
start out with a similar slope as in Fig. 



terms of the lattice interaction U = g c gn. The tran- 
sition curve is now determined by the implicit equa- 
tion Fjj (ksT / J,U I J) = 0, where in zeroth hopping or- 
der F^ 0) (x, y)=x - 2y/D 2 + Dy/ ln^D 2 + Dy + y + 
2D]/[A^D 2 + Dy - y - 2D]. The resulting transition 
curve for D — 3 and the next four approximations com- 
ing from successive hopping orders are shown in Fig. [21 
A fast convergence is observed, with the approximation 
sequence of transition points at T — corresponding to 
(U/J)T=° = 6(3 + 2%/3) w 38.8, 34.1, 32.2, 31.8, 31.6 . . . 
which converge to roughly 30.8, as shown in Fig. [3] 
Thus our value is slightly smaller than the mean-held rc- 
su lt (U/J) ? =0 ~ 34.8 derived from Bose-Hubbard model 
[H l3ft l3ll l32j | and the experimental number (C//J)^ =0 « 
36 |5j- The associated hopping sequence of transition 
temperatures at U — converges to w 3.6 J/ks- 
For the lattice spectrum (1) we cannot improve the 

result near T c (0) in the same way as for the free-particle 
spectrum. By analogy, we may, however, assume that the 
characteristic reentrant transition will also here survive 
higher-loop corrections. 

An important question is whether the nose survives the 
experimental set-up of an optical lattice where the par- 
ticles require an extra weak overall magnetic trap to 
keep them together. Indeed, if the magnetic trap is 
harmonic as in all present experiments with typical fre- 
quency Wtrap ~ 2ir x 24 Hz, the nose is expected to 
disappear since by analogy with a continuum calcula- 
tion, the external trap reverses the slope of the transition 
curve at tJ - 1 [33l l34l | (see also Chapter 7 in the textbook 
However, a nose will definitely be seen if the exper- 
iment is performed with a magnetic trap which has an 
approximate box-like magnetic shape. This can be ex- 
perimentally realized by using optical dipolar traps and 
lasers with square profiles. It will be interesting to study 
experimentally the dependence of the slope of the tran- 
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FIG. 3: Convergence of hopping expansion for the critical 
value of U I J at the zero-temperature quantum phase transi- 
tion. The limiting value for (U/J)^ =0 is 30.8. 



sition line on the shape of the magnetic trap potential as 
it is distorted from harmonic to a box-like shape. For a 
power-like shape of the trap |85l th ere should be a critical 
power where the nose appears |36|. 
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